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Abstract 

We interpret magnetic billiards as Finsler ones and describe an 
analog of the string construction for magnetic billiards. Finsler bil- 
liards for which the law "angle of incidence equals angle of reflection" 
are described. We characterize the Finsler metrics in the plane whose 
geodesies are circles of a fixed radius. This is a magnetic analog of 
Hubert's fourth problem asking to describe the Finsler metrics whose 
geodesies are straight lines. 
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1 Introduction and background material 

This paper concerns the motion of a charged particle in a magnetic field, a 
popular object of study in mathematics and mathematical physics. In the 
Euclidean plane, the strength of the magnetic field is given by a function 
B{xi,X2), and the particle moves with a constant speed, satisfying the equa- 
tion 

x = B{xi,X2)Jx where J(f i, ^2) = (— f2, f i) (1) 
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If one fixes the speed \v\ then the magnetic field prescribes the curvature of 
the trajectory at every point. In particular, if the field is constant then the 
trajectories are circles of the Larmor radius |f Our sign convention is 

that if -B > then the circles are traversed in the counterclockwise direction. 

In general, a magnetic field on a Riemannian manifold M is a closed 
differential 2-form j3, and the magnetic fiow is the Hamiltonian fiow of the 
Riemannian Hamiltonian function |pp/2 on the cotangent bundle T*M with 
respect to the twisted symplectic structure u + Tr*{P) where u = dpAdq is the 
standard symplectic structure on T*M and vr : T*M — >■ M is the projection. 
We refer to jH [T31 US UHl CHI for a variety of results on magnetic flows on 
Riemannian manifolds. 

If a charged particle is conflned to a domain with ideally reflecting bound- 
ary then one has a magnetic billiard. The particle moves inside according to 
equation ((T)) and undergoes elastic reflections off the boundary: the tangen- 
tial component of the velocity remains the same and the normal one changes 
sign. In dimension two, this amounts to the familiar law of geometrical op- 
tics: the angle of incidence equals that of reflection. Magnetic billiards have 
attracted a considerable attention: see [TUl HH [121 UHl I2ni 1211 12H1 IBS] ; see also 
for a survey of various aspects of billiard systems. 

In this paper we interpret a magnetic flow as a geodesic flow of a Finsler 
metric. We mostly consider the 2-dimensional case. In Section 2 we interpret 
magnetic billiards as Finsler ones and describe the magnetic version of the 
string construction that recovers a billiard table by a caustic of the billiard 
map. We also characterize the Finsler metrics for which the Finsler billiard 
enjoys the familiar law "angle of incidence equals angle of reflection" (The- 
orem |21 and Corollary El) . In Section 3 we describe the Finsler metrics in 
the plane whose geodesies are circles of a flxed radius; we give analytic and 
synthetic descriptions in Theorems I6I7I81 This is an analog of the celebrated 
Hubert's fourth problem of describing the Finsler metrics in a domain in pro- 
jective space whose geodesies are straight fines (H lU (20]. Our solution 
has an unexpected connection to another classical and well studied question: 
the Pompeiu problem [301 1^ IH^. 

We will now review basics of Finsler geometry (see, e.g., j21 El El I21j) and 
Finsler billiards, recently introduced in ^7]. Finsler geometry describes the 
propagation of light in an inhomogeneous anisotropic medium. This means 
that the velocity of light depends on the point and the direction. There are 
two equivalent descriptions of this process corresponding to the Lagrangian 
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and the Hamiltonian approaches in classical mechanics, and we will mostly 
use the former. 

The optical properties of a medium are described by a quadratically con- 
vex smooth hypersurface, called the indicatrix, in the tangent space at each 
point. The indicatrix consists of the velocity vectors of the propagation of 
light at a point in all directions. It plays the role of the unit sphere in Rie- 
mannian geometry. We do not assume that the indicatrices are centrally 
symmetric. 

Equivalently, a Finsler metric on a manifold M is determined by a smooth 
nonnegative fiber-wise convex Lagrangian function L[x,v) on the tangent 
bundle TM, homogeneous of degree 1 in the velocity: 



The restriction of L to a tangent space T^M gives the Finsler length of 
vectors in T^M, and the indicatrix at x is the unit level hypersurface of 
L{x, v). Given a smooth curve 7 : [a, b] — >• M, its length is 



The integral does not depend on the parameterization. A Finsler geodesic is 
an extremal of the functional C The Finsler geodesic flow is the flow in TM 
in which the foot point of a vector moves along the Finsler geodesic tangent 
to it, so that the vector remains tangent to this geodesic and preserves its 
Finsler length. The Finsler geodesic flow is described by the Euler-Lagrange 
equation 

dLy{x,v)/dt — Lx{x,v) = or L^yV + Ly^v — = 0. (2) 

The dual, Hamiltonian approach describes the propagation of light in 
terms of wave fronts and the Finsler geodesic flow as a Hamiltonian flow in 
the cotangent bundle T*M. Let / C T^M be the indicatrix. The figuratrix 
J C T*M is the dual hypersurface constructed as follows. Given a vector 
u & I, the respective covector p e J is defined by the conditions: 



This gives a diffeomorphism I ^ J, called the Legendre transform. In the 
same way as the field of indicatrices determines the Lagrangian L, the field of 



L{x,tv) = tL{x,v) for all t > 0. 




Ker p = T^I and p{u) = 1. 
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figuratrices determines a Hamiltonian H on T*M. The Hamiltonian vector 
field of tlie function H is also called the Finsler geodesic flow; the Legendre 
transform identifies the two flows. 

Example 1: Hilbert's fourth problem in dimension two. The Eu- 
clidean metric is given by the Lagrangian L{x,v) = \v\] its geodesies are 
straight lines. Such metrics are called projective. Following [T], let us de- 
scribe all symmetric projective Finsler metrics in the plane, that is, a solution 
to Hilbert's fourth problem in dimension 2. 

A synthetic approach, due to Busemann, makes use of integral geometry, 
namely, the Crofton formula j22]- Consider the set of oriented lines in the 
plane, topologically, the cylinder. An oriented line can be characterized by 
its direction a G [0, 27r) and its signed distance p from the origin. The 2-form 
ujq = dp A da is the standard area form on the space of oriented lines; this 
symplectic form is a particular case of a symplectic structure on the space of 
trajectories of a Hamiltonian system on a fixed energy level, in particular, the 
space of oriented geodesies of a Finsler metric - see, e.g., [71 and a discussion 
in Section 3. The Crofton formula gives the Euclidean length of a plane curve 
7 in terms of Uq. The curve determines a function on the space of oriented 
lines, the number of intersections of a line / with 7. Then 

length(7) = (1/4) J #(/ n 7) cuq. (3) 

Let f{p,(y) be a positive continuous function. Then u = f{p,a) dp A da 
is also an area form on the space of oriented lines. Formula 0, with lj 
replacing c^o, defines a projective Finsler metric, and all such metrics can be 
obtained by an appropriate choice of the function /. 

Next we describe an analytic solution to Hilbert's fourth problem in di- 
mension two. First, a Lagrangian L{x, v), homogeneous of degree 1 in v, gives 
a projective Finsler metric if and only if the mixed second partial derivative 
matrix Lxv is symmetric; this is Hamel's theorem of 1903, and it holds in any 
dimension. The Lagrangians satisfying Hamel's condition have the following 
integral representation: 

L(xi, X2, f 1, ^2) = / |fi COS0 + f2 sin0| /(xi COS0 + X2 sin0, 0) (i0 (4) 
Jo 

where f{p,(f)) is a smooth positive function on the cylinder representing the 
space of oriented lines. Moreover, if / is even in (p then it is uniquely deter- 
mined by L. The function / is the same as in Q: the length of a curve with 
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respect to the Finsler metric is given, up to a multiplicative constant, by 
0. If / depends on the angle only then one obtains a translation invari- 
ant metric, called a Minkowski metric. If / is a constant then one has the 
Euclidean metric. 

Let M be a 2-dimensional Finsler manifold with boundary, a curve N. 
The Finsler billiard system is defined in as follows. A point moves inside 
M freely, according to the Finsler geodesic flow, until it hits the boundary. 
The reflection is described in terms of the indicatrix I at the impact point 
X - see figure 1. The vectors u and v are the Finsler unit vectors along 
the incoming and outgoing trajectories. The tangent lines to / at u and 
V are concurrent with the tangent line to at x. This definition satisfies 
a variational principle: for every points a,b & M, the reflection point x 
extremizes the Finsler length |ax| + If the indicatrix is a circle centered 
at the origin then the vectors u and v make equal angles with the boundary 
curve N; this is the familiar law of Euclidean billiard reflection. The multi- 
dimensional version of the Finsler billiard reflection is defined similarly, and 
we do not dwell on it - see [T7j. 




Figure 1: Finsler billiard reflection 

Example 2: projective Finsler billiard reflection. Consider a sym- 
metric projective Finsler metric (^. In polar coordinates, vi = rcosa,V2 = 
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rsina, and 

/■Q!+7r/2 

X2, r, a)/2 = r / cos(a — 0) /(xi cos0 + X2 sin0, 0) (5) 

J a--K/2 

Let a be the direction of the bilhard curve at the impact point x and (3 and 
7 the directions of the incoming and the outgoing bilhard trajectories - see 
figure 2. The projective Finsler reflection law specializes to the following 
formula. 

Lemma 1.1 One has: 

/ cos(a — (f)) f{xi cos (j) + X2 sin 0, 0) dcf) = 

J'y-TT/2 
/3+7r/2 

cos(q; — 0) /(xi cos + X2 sin 0, 0) (i0. 

+7r/2 




Figure 2: Deriving the projective Finsler billiard reflection law 

For example, if / = 1 then integration yields: cos(/3 — a) = cos(7 — a) or 
/3 — a = 7 — a, the familiar law of equal angles. 

Proof. Denote the integral in by g{x, a). Then the polar equation of 
the indicatrix at point x, chosen as the origin, is r = 1/ g{x, a). It is a matter 
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of a straightforward calculation to find the coordinates of the intersection 
point P in figure 2: 

^^^^ (cos a, sin a) 



g{(3) cos(a - (3) + g'{(3) sm{a - f3) ' 
Equating P{P) and -P(7) yields: 

g{P) cos(a - (3)+ g'{/3) sin(a - /3) = gin) cos(a - 7) + ^''(7) sin(a - 7). (6) 
It follows from (jH)) that 

'/3-7r/2 

and similarly for 7. It remains to substitute to ^ and to collect terms. □ 



g'{l3) = — sin(/5 — 0) /(xi cos + X2 sin i 

J/3-7r/2 



2 Magnetic billiards as Finsler billiards 

Consider the plane motion of a charged particle in a magnetic field with 
strength B{xi,X2)- The Lagrangian for this motion is 

L{x,v) = ^\v\'^ + a{x){v) 

where a{x) = f{xi,X2) dxi + g{xi,X2) dx2 is a differential 1-form such that 
da = —B{xi,X2) dxi A dx2- The choice of a is not unique: one can always 
add a closed 1-form to a Lagrangian without effecting the dynamics. The 
Euler-Lagrange equation for L is (P). In particular, the Lagrangian for a 
constant magnetic field is 

L{x,v) = + ^[v,x] 

where [ , ] is the cross-product. 

Following the Maupertuis principle (see, e.g., 0), we replace the La- 
grangian L by 

L{x,v) = \v\ + a{x){v). (7) 
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The extremals of the Lagrangian (|7j) coincide with those of L, corresponding 
to the motion with the unit speed. In particular, the extremals of 



v) 




(8) 



are the counterclockwise oriented circles of radius R. 

The Lagrangian ((Tj) defines a non-symmetric Finsler metric in the domain 
where L{x,v) > for all f 7^ 0. This is the case if \a{x)\ < 1, and we 
assume this condition to hold throughout this section. In other words, we 
assume that the magnetic field is sufficiently weak. Under this assumption, 
we consider the unit speed magnetic flow as the Finsler geodesic flow. 

Consider a plane domain and the magnetic billiard inside it. One also has 
the Finsler billiard inside the domain, associated with the Lagrangian ((7j). 
One expects the two systems to coincide, that is, to have the same reflection 
laws. 

Theorem 1 The Finsler billiard reflection law, associated with the Lagrangian 
is the law of equal angles: the angle of incidence equals the angle of re- 



Proof. The indicatrix of the Finsler metric at point x is given by the 
equation |t>| +a{x){v) = 1. Choose Cartesian coordinates in such a way that 
a{x){v) = tvi for some t G R. Then the equation of the indicatrix can be 
rewritten as 



Recall that the equation of a conic, centered at the origin in the {vi, f2)-plane, 
is 



where — = h"^. Clearly, Q has this form. Hence the indicatrix is an 
ellipse, centered at the origin. 

Thus the theorem reduces to the following geometrical property of conies. 
Let / be a conic with focus O, let X, F G /, and let Z be the intersection 
point of the tangent lines to / at X and Y . Then the line OZ bisects the 
angle XOY - see figure 3. This property holds indeed; it is known as the 
Poncelet "first little theorem", see This completes the proof. □ 



flection. 




1. 



(9) 
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X 




In fact, the law of equal angles is characteristic of the conies, centered at 
a focus. 

Theorem 2 Let I be a smooth plane curve, star-shaped with respect to point 
O, with the following property. Let X,Y & I be arbitrary points, Z be the 
intersection point of the tangent lines to I at X and Y ; then the line OZ 
bisects the angle XOY. It follows that I is a conic with focus O. 

Proof. Let O be the origin, and give / a parameterization I{t) so that 
[I{t), r{t)] = 1. Then J"(t) = -f{t)I{t) for some function f{t); thus we view 
/ as an orbit in a central force field. We claim that f{t) = C/\I{t)\^ for a 
constant C. Assuming this claim, it follows that I is an orbit in Newton's 
force field, and therefore a conic with focus O. 

Let X = I{ti) = Ii,Y = 1(^2) = h- A direct computation yields the 
point Z: 

The equal angle condition reads: \Y,Z]/\Y\ = [Z, X]/\X\, or 

(l-[/i,J^]) = |/i| ii-[h,i[]). (10) 

Now set: ti = t,t2 = t + e and shorthand I{t) to I and f{t) to f. Then the 
Taylor expansion yields: 
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and 

I/2I = |/|+£^ + 0(£2). 

Substitute to (fTUI) and collect terms to obtain: + 3/ ■ /'/ = 0. This 

differential equation is easily solved: /'// = —3/ ■ /'/l-^P hence / = 
C/|/p, as claimed. □ 

As a consequence, we obtain a description of Finsler metrics for which 
the Finsler billiard reflection law is the law of equal angles. 

Corollary 3 The Finsler billiard reflection satisfies the law "angle of in- 
cidence equals angle of reflection" for every billiard curve if and only if the 
metric is given by a Lagrangian 

L{x,v) = f{x){\v\ + a{x){v)) (11) 

where f{x) is a non-vanishing function and a{x) is a 1-form. 

Proof. Replacing a metric by a conformally-equivalent one changes the 
indicatrices by a dilation and does not effect the law of equal angles. Theorem 
^implies that the metrics (fTT|) satisfy this law of equal angles. Conversely, 
if this law holds then, by Theorem |21 the indicatrices are ellipses (depending 
on the point of the plane), centered at their foci. A general equation of such 
an ellipse is f{x){\v\ + a{x){v)) = 1, and the result follows. □ 

It is shown in jT7] that some familiar properties of the usual billiards ex- 
tend to the Finsler ones. Although |2| concerned symmetric Finsler metrics, 
the results hold in the non-symmetric well; however one should be 

careful with the order of points: the distance from A to B may differ from 
the distance from B to A. Let us consider the case of a constant magnetic 
field, that is, the Finsler metric given by the Lagrangian (jSJ whose geodesies 
are counterclockwise oriented arcs of radius R. 

Let A and B be two points on an arc of radius R with the center C and 
the angle measure 6. Denote the Finsler distance between points by d{A, B) 
and identify points with their position vectors. Let >C(7) denote the Finsler 
length of a curve 7. 
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Lemma 2.1 One has: 

d{A,B) = h)R + ^[B-A,C]. 
For a simple oriented closed curve 7, one has: 

^(7) = /(7) - ^5(7) (12) 

where l{'~f) and S{'j) are the Euclidean length and the Euclidean signed area 
bounded by 7. 

Proof. To obtain ()12j) . one integrates |f | + [v,x]/{2R) over 7 and makes 
use of the fact that [v, x]/2 is negative the derivative of the signed area swept 
by the position vector of 7. 

Let O be the origin. The distance d{A,B) equals the integral of |f| + 
[v,x]/{2R) over the arc AB. The integral of |f| is the arclength of the arc, 
that is, OR. The integral of [v,x]/{2R) equals —S/R where S is the area of 
the curvilinear triangle OAB. 

Assume first that O = C. Then the latter area is 6'i?^/2, and d{A.,B) = 
6R/2. If the origin is translated through vector C then the Lagrangian 
changes by the term [v,C]/{2R), and its integral by [B — A,C]/{2R). This 
yields the first formula. □ 

Remark 2.2 Formula ()12j). along with its proof, holds for closed immersed 
curves as well: the area term should be understood as the integral of the 
1-form {xdy — ydx) /2 over the curve. 

The orientation of 7 determines a coorientation: the pair (coorientation 
vector, orientation vector) gives the positive orientation of the plane. If 7 
is a counterclockwise oriented simple curve then the positive coorientation 
is the outward one. Given a real number t, consider the parallel curve r{t) 
at distance t from 7. The curve r{t) is the time-t wave front, starting at 
7 = r(0). More precisely, one translates the contact elements of 7 in the 
orthogonal direction through distance t (along the coorienting vector, if t > 0, 
and in the opposite direction, if t < 0), and the obtained 1-parameter family 
of contact elements consists of the contact elements of T{t). The curve r[t) 
may have singularities, generically, semi-cubic cusps. If 7 is a positively 
oriented circle of radius r then r{t) is a circle of radius r + t. 

Formula (jl2p admits the following interpretation. 



11 



Lemma 2.3 For an oriented simple closed curve 7, one has: 

C{^) = ^{nR'-S{T{-R))) (13) 

where the area S{r{—R)) is understood as in B.emark \2.^ 

Proof. For a closed immersed curve, the following well known formula 
holds: 

SiTit)) = Si^)+tli^) + nt^w 

where w is the Whitney winding number of 7. Therefore the right hand side 
of (fT^ equals (l/-R)(7r_R^ — S{T{—R))), and the result follows from Lemma 

o □ 

Remark 2.4 Formula (jl3|) expresses Finsler lengths in terms of areas, and 
it serves a magnetic analog of the Crofton formula (jS)). The curve r{—R) is 
the locus of the centers of positively oriented circles of radius R, tangent to 
7 and having the same orientation as 7 at the tangency point. The curve 7 
can be reconstructed from r{—R) as the envelope of the family of circles of 
radius R centered at points of T{—R). 

Formula ()13j) also resembles the area-length duality for spherical curves, 
discussed in 1211 121] • Let 7 be a simple smooth closed curve on the unit 
sphere, and let F be its spherically dual curve, namely, the curve T{n/2), in 
the sense of spherical geometry. Then one has: /(7) = 27t — S{r). 

Let us return to billiards. Recall that a caustic of a 2-dimensional billiard 
is a curve F inside it with the following property: if a segment of a billiard 
trajectory is tangent to F then so is the reflected segment. Given a convex 
caustic F, can one reconstruct the billiard table? For the usual, Euclidean 
billiard the answer is given by the string construction: a billiard curve N is 
the locus of points X of a string of fixed length, wrapped around F - jHl 123 • 

It is shown in ^7] that the string construction extends to Finsler billiards 
as well. For non-symmetric Finsler metrics one needs to consider oriented 
caustics F so that the orientation of billiard trajectories, tangent to F, agrees 
with the orientation of F. Applying these considerations to billiards in a 
constant magnetic field, we obtain the following corollary. 

Let F be an oriented closed convex curve. For a point X outside of F, let 
F{X) be the Finsler length, associated with the metric (jH)), of the shortest 
closed curve from X to X around the obstacle F, whose orientation agrees 
with that of F - see figure 4. 
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X 




Figure 4: Magnetic string construction 



Corollary 4 The level curves of the function F{X) are the boundaries of 
magnetic billiard tables that have T as a caustic. 

Recall the optical property of an ellipse: a ray emanating from one fo- 
cus reflects to another focus. As a particular case of Corollary El one may 
construct a magnetic analog of an ellipse. 

Corollary 5 Let A and B be fixed points ("foci") and N be the locus of 
points X such that d{A, X) + d{X, B) = const. Then every trajectory of the 
magnetic billiard, starting at A, reflects in N to B. 

Note that the two foci play different roles and cannot be interchanged in 
the above formulation. If the points A and B merge then the "ellipse" 
becomes a Euclidean circle centered at this point. 

3 Finsler metrics whose geodesies are circles 
of a fixed radius 

In this section we develop a magnetic analog of the solution to Hilbert's fourth 
problem, outlined in Example 1 above. Start with an analytic description of 
the Lagrangians, homogeneous of degree 1 in the velocity, whose extremals 
are positively oriented circles of radius R. 
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Lemma 3.1 The extremals ofL{x,v) are positively oriented circles of radius 
R if and only if L satisfies the equation: 



1^1 

'r 

where J(fi,f2) = (— f2,fi 



Lvv{Jv) + L^x{v) = (14) 



Proof. Let x{t) be a parameterized curve, v = x' . The curve is a 
counterclockwise oriented circle if and only if 



V \' 1 



Differentiate to express the acceleration vector: 

^ = y ^(^) + (15) 

Since the Lagrangian is homogeneous of degree 1, the Euler equation L^v = L 
holds, and hence L^viv) = 0. It remains to substitute v' from (fTH|) to the 
Euler-Lagrange equation 

and the result follows. □ 

We are ready to prove the main analytical result of this section. 

Theorem 6 Every Lagrangian, homogeneous of degree 1 in the velocity, 
whose extremals are positively oriented circles of radius R can he represented, 
in polar coordinates, as follows: 

L{x,v) = L{xi,X2,r,a) = ri / cos{a—(j)) g{xi+Rcos(f), X2+Rsm(f)) d(f) + 

a{xi,X2) cosa + b{xi,X2) sina^ (16) 

where g is a positive density function in the plane such that the center of 
mass of every circle of radius R is its center, and a, b are two functions, 
satisfying 

(3:^1, 3:2) - bxi{xi,X2) = 4 9{xi + R,X2). (17) 
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Proof. In polar coordinates, Vi = rcosa, V2 = rsina, and one has: 
L{x,v) = \v\ p{xi,X2, a) for some function p. 

Fix a point x = {xi,X2) and consider the indicatrix / at x, chosen as the 
origin. The polar equation of / is r = l/p{x,a). Therefore p{x,a) is the 
support function of the dual curve, the figuratrix J (see |22j). Parameterize 
J by the angle made by its tangent vector with the horizontal axis. Let 
f{x,(l)) be the radius of curvature at point J(0) and let {a{x),b{x)) be the 
coordinates of the point J(0). One has: 

J'i(l>) = /(2;,0)(cos0,sin0), 

and hence 

/■a+7r/2 

J(a + 7r/2) = J(0) + / /(x,0)(cos0,sin0) c?0, 

Jo 

see figure 5. It follows that 

p{x, a) = (cos a, sin a) ■ J{a + n/2) = 



/•a+TT/2 

a{x) COS a + b{x) sin a + / cos{a — (f))f{x, (f)) dcj). (18) 

Jo 

Differentiating (|18|) twice, one recovers the function / from p: 

f{x, a + n/2) = p{x, a) + p"{x, a). 

Note that every function of the form p + p" is orthogonal to cos a and 
sin a. Thus 

/'2n p2n 

/ f{x, a) cosa da = / f{x, a) sin a da = 0; (19) 
Jo Jo 

this also follows from the integral representation p8|) and periodicity of p as 
a function of a. 

Next we use the equation (fT^ in the integral representation (fTHj) . One 
rewrites the differential operators d'^v, dvdx and dx in polar coordinates and 
applies to L{xi,X2,r,a) = rp{xi,X2,a), given by (fTH|) . Taking into account 
that V = r(cosQ;, sina), a computation reveals that 

\v\ rf{x,a + 7i/2) . 

—K.,[.Jv) = (-sma,cosa), 

rC ri 
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Figure 5: Integral representation of the support function 



and 



- L^^{v) = r{ I (cos {x,(f)) - sin (p f^^ {x,(f))) 



+ 



Therefore, by ()14j] 
1 



0'X2ix) — bx-^{x) ] (— sin a, cos a) 



Q+7r/2 



^f{x,a + Tr/2) = j (cos0 fx^^x^cj)) - sin0 /^^(x, 



+ 



In particular, 



1 /(x, 0) = a^2(x) - 



Differentiating ipUl) with respect to a, one gets: 

4 faix,a) = cos a fx2{x,a) - sin a fxj_{x,a). 
K 
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We claim that 



/(xi, X2, a) = g{xi + R cos a,X2 + R sin a) (23) 

for an appropriate function of two variables g. 
Indeed, consider the vector field 

T] = — da + sin a dxi — cos a 8x2 
R 

on the solid torus x S^. Then (j22|) can be written as rj{f) = 0. The 
trajectories of 77 are: 

a{t) = xi{t) = -i?cos (^-^ + ^2(i) = --Rsin ^-^ 

where t is the "time" parameter and a, b are constants. One can take the 
plane a = as a section. Then the T^-invariant function / is determined by 
its values on this section, a function of two variables g. Consider a point 
{a,Xi,X2). The trajectory through this point intersects the section at point 
(a, 6) = (xi + -Rcosa;,X2 + -Rsina). Hence /(xi,X2,a) = g{a,b), and 
follows. 

Equations (fTU|) imply that the center of mass of the circle of radius R, with 
the density function g, centered at {xi,X2), is the point {xi,X2). Combining 
(EH) and (Uni), we obtain (HH). □ 



Remark 3.2 The term r{a{x) cos a + b{x) sin a) in ()16|) . the formulation of 
Theoremini can be written as i'{x){v) where z/(a;) = a{x) dxi+b{x) dx2 is a 1- 
form. The choice of this form is not unique but dv is uniquely determined by 
the function g via (jl7|) . This is consistent with the remark we already made: 
adding a closed 1-form to the Lagrangian does not effect the Euler-Lagrange 
equations. 

To proceed, we recall basic facts about the symplectic reduction. Let 
(M, VL) be a symplectic manifold and if : M — > R a Hamiltonian function. 
Consider the Hamiltonian vector field ^ = sgrad H . Since H is ^-invariant, 
the field ^ is tangent to the level hypersurfaces of H. Consider such a hyper- 
surface S", and assume that the space of trajectories of ^ on is a smooth 
manifold = S/^\ locally, this is always the case. The restriction of to 
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S has a 1- dimensional kernel spanned by ^, and hence descends to a 
symplectic structure u on A^. This is the symplectic reduction of Q. 

One applies this construction as follows. Given a Finsler manifold M, 
the symplectic manifold in question is the cotangent bundle T*M with its 
standard symplectic structure dp A dx where x G M is the position and 
p &T*M the momentum. The function H is the Finsler metric Hamiltonian, 
and the hypersurface S consists of the unit covectors; it is fibered over M and 
the fibers are the figuratrices. The vector field ^ is the Finsler geodesic flow, 
and the space of trajectories identifies with the space of non-parameterized 
oriented geodesies. 

Consider the tangent bundle TM and the unit vector hypersurface U in 
it. The Legendre transform (x, f) ^— = L„) identifies U with S and the 
Finsler geodesic flow on U with the geodesic flow ^ on S. The pull-back 
of the Liouville form pdx is the 1-form A = L^dx on TM. The form A is a 
contact form on f/, and C is its Reeb vector field: A(C) = 1, ic^dX = 0. The 
reduction of the 2-form dX yields the symplectic structure on the quotient 
space U/(, the space of oriented Finsler geodesies. 

Given a smooth curve 7 on M, one lifts it to the curve 7 on [/ by assigning 
the unit tangent vector to every point of 7. Then the Finsler length of 7 
equals 



For a reference on this symplectic approach, see, e.g., jSlEl- 

Now we are in a position to compute the symplectic structure on the 
space of circles of radius R, associated with the Lagrangian (fT^ . A circle is 
characterized by its center, and the space of circles is the plane with Cartesian 
coordinates {u, v). 

Theorem 7 The symplectic structure uo on the space of circles of radius R, 
associated with the Lagrangian 17^) . is given by the formula: 



Proof. The manifold U consists of the Finsler unit tangent vectors in 
the plane and has coordinates a,xi,X2. We use the notation from the proof 
of Theorem IHl The formulas derived in that proof yield: 




g{u, v) du A dv. 



(24) 
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f a+7r/2 

A = Lydx = I / COS0 f{x,(l)) dcj) + a{x) 1 dxi + 



'0 



' /.a+7r/2 \ 

/ sincj) f{x,(j)) d(j) + b{x) 1 dx2 



and, taking and (PT|) into account, 

= /(x, a + 7r/2) f cosa lia A (ix2 — sin a da A dxi — '^^xi A (ia;2 ) . 



In view of 

dX = g{xi — R sin a,X2 + R cos a) ^cos a rfa A (ix2— 

sin a da A dxi — -dxi A 6^X2) . (25) 
R ' 

Now consider the projection U U / Q = R^. To compute the symplectic 
structure in R^, consider a section j : R^ — » t/ and let uj be the pull-back 
of d\] the result is independent of the choice of j. As a section one may take 

TT 

j{u, v) = {a, xi, X2) with ^ — 1^^ xi = u + R, X2 = v, 
see figure 6. It remains to substitute to (j^ . and the result follows. □ 

As a consequence, the Finsler metric (fT^ can be recovered, up to sum- 
mation with a closed 1-form, from the area form (|24j) - see Remark 13.21 

Next we consider an analog of formula (jl3p for a general Finsler metric 
(fTB|) whose geodesies are circles of radius R. The following result expresses 
the Finsler length in terms of the area form on the space of circles and is 
analogous to the synthetic solution to Hilbert's fourth problem, that is, the 
Crofton formula, used as a definition of a projective metric. 

Theorem 8 Given an oriented simple closed curve 7, one has: 

C{^) = S{T{-R)) + C (26) 

where S{r{~R)) is the area bounded by the curve r(— i?) with respect to the 
area form ^24\ ), and C is the common Finsler length of all positively oriented 
circles of radius R. 
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(u,v) 




Figure 6: Computing the area form on the space of circles 



Proof. Note first that the geodesies are extremals of the length functional 
L. The space of geodesies identifies with the plane which is a critical manifold 
of L. A function is constant on its critical manifold, hence all positively 
oriented circles of radius R have equal Finsler length. 

To prove the result we consider a variation of the curve 7 and show that 
both sides of (j26|) have the same variations. This being established, one can 
deform 7 to a circle of radius R for which the result holds. 

Assume that 7 is parameterized by the Euclidean arc-length. Then 7' = 
(cos a, sin a). The lift of 7 to t/ is the curve 7 = (7, a). Consider a variation 
of the curve, that is, a vector field w along 7. It is straightforward to compute 
that the respective variation of 7 is the vector field 



where A is the contact form as in the proof of Theorem [7| Therefore the 
variation of the length £(7) is given by the formula 



w = w \ [7', v'\ da. 



(27) 



One has: 




7 




7 
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where dX is as in (j25|) . Let k{t) be the curvature at 7(t). Then da = kdt. A 
computation using and (|77j) reveals that 

Jiu,dX = J g{xi — Rsma,X2 + Rcosa) [j^w] (^'^^^^^- (^^) 

On the other hand, one has: 

r(— i?) = (Xi, X2) = (xi — _Rsina, X2 + Rcosa). 
Therefore the variation of r(— i?) is given by the vector field 
u = {wi — Riy, v'] cos a, W2 — Riy, v'] sin a). 

Then 

dXi = (1 — Rk) cos a dt, dX2 = (1 — Rk) sin a dt. 

Since u = (— l/i?) g{Xi,X2) dXiAdX2, it is straightforward to compute the 
variation of the area S{T{—R)): 

f iu^ = 1 g{xi — Rsina, X2 + Rcos a) [y\w] ( — k] dt. 
Jr{-R) J \^ J 

This is the same as (EHI), and we are done. □ 



Note the following corollary of formula ()26p . 



Corollary 9 The integral of the area form \24^ is the same over all discs of 
radius R. 

Proof. Let 7 degenerate to a point in ()2(ij) . so that C{'~f) = 0. Then 
T{—R) is a circle of radius R, and the cu-area, bounded by it, equals —C. □ 



Remark 3.3 One can give a somewhat different proof of Theorem El that 
does not use the specifics of the Euclidean plane and applies to other surfaces, 
for example, the sphere. Let us outline the argument. Pick a point O inside 
7 and consider an infinitesimally small loop S around O whose orientation 
is the same as that of 7. The Finsler unit tangent vector fields to 7 and to 
6 extend to a unit vector field in the annulus A bounded by 7 and 6. This 
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vector field provides a lift A of the annulus to U, and dA = 'j — S. By Stokes' 
theorem, 



A - / A = d\. 

The second integral on the left hand side is infinitesimally small. The integral 
on the right can be understood as the symplectic area of the set of circles of 
radius R whose centers he between the curves T{—R) and the circle of radius 
R, centered at O, and ()2fi|l follows. This also shows that the symplectic area 
of a circle of radius R is independent on its choice. 

One can revert the arguments and and define the respective Finsler met- 
ric, as in Theorem ini starting with an area form u = g{xi,X2) dxi A dx2i 
satisfying the property that the w-area of every disc of radius R is the same. 
Then the function g{xi^X2) should be orthogonal to cosine and sine on 
every circle of radius R. 

Lemma 3.4 The integrals of a function g over all discs of radius R is the 
same if and only if g is orthogonal to cosine and sine on every circle of radius 
R. 

Proof. Let S* be a circle of radius R with center x = {xi,X2)- Consider 
its variation given by an infinitesimal parallel translation through vector 
V = (wi, ^2). The variation of the cj-area of the disc is 

/ iyUJ = / g{xi + R COS a, X2 + Rsin a) {vi sin a — V2 COS a) da. 
Js Jo 

This vanishes for all v if and only if g is orthogonal to cos a and sin a. □ 

How restrictive are these two equivalent conditions on function g7 This 
question goes to the heart of the Pompeiu problem, see jSOJ ED 112] • Given 
a compact set K, one considers the continuous functions with zero integrals 
over all isometric images of K. For which sets K must such functions be 
identically zero? D. Pompeiu, who posed this problem in the late 1920-s, 
erroneously thought that the disc in the plane has this property. In fact, 
if i^' is a disc then there are plenty of functions with zero integrals over all 
congruent discs; although there is a wealth of results on this subject, the 
general solution to the Pompeiu problem is not known yet. 

If a function g has a constant integral over all discs of radius R then it 
can be written as Const + h{xi,X2) where h has zero integrals over all discs 
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of radius R. The following result, standard in the literature on the Pompeiu 
problem, provides a substantial supply of such functions. We need to recall 
the definition of the Bessel functions. 

The Bessel functions J„(tf ), n G Z, are defined by the generating function 



\ \ / / n— — rv^ 



(29) 



An exphcit formula is as follows: 

■^"H = Z^ 2"+2j(n + j)!j! ' -^-"H = ■^"(-^)- 

We will need the following property: 

y Jo(w) wdw — wJ\{w). (30) 

Lemma 3.5 Let a he a root of the first Bessel function Ji, and let f{(3) be 
a function on the circle. Then the functions 

h{xi,X2) — J cos ^-^ (a;iCOS/3 + a;2sin/3)^ f{f3) df3 (31) 

and ^ 

h{xi,X2) = J sin ^-^ (xiC0S/3 + X2sin/3)^ f{P) dp 

have zero integrals over all discs of radius R. 

One may also take linear combinations of such function over different 
roots of Ji. 

Proof. Let D be the disc of radius R centered at the origin, and let ^ 
be its characteristic function. The condition on function h reads: ^ * /i = 
where * denote the convolution. Take the Fourier transform to obtain: 

I /i = 0. (32) 

Let us compute ^: 

^^(A) = / e"'^" dx= / e-''''"'°<''-^^ da rdr 

J D Jo Jo 

23 



where x = r (cos a, sin a), A = p(cos/3, sin/3). One has: 

-trpcos{a- (3) = ^ {e'' - e-'') 

where 9 = a — /3 — 7r/2. Using the definition of Bessel functions ^I^. it 
follows that 

/ e-^"^™'^^"-^) da = y] Jn{rp) / e*"^ dO = 27rJo(rp). 
Jo „ Jo 

By (jHUjl . one has: 

/ Jo{rp)rdr = —Ji{Rp) 
Jo P 



and hence 



e(A) = ^Ji(i?|A|). 



The condition ()H2|1 holds if the support of h is contained in the union of 
circles, centered at the origin, whose radii are a/R where a is a root of Ji. 
Fix one such root and let 

= /(/?) 5.(p). 



Taking the inverse Fourier transform yields: 

p27r POO 



Jo 



a 
R 



2tt 

One concludes by taking the real and imaginary parts. □ 



For example, let / in (p?T|l be the delta function 6o- Then h = cos{axi/ R). 
Substitute g = 1 + cos{axi/ R) into to obtain an "exotic" Finsler metric 
whose geodesies are circles of radius R. 

Remark 3.6 One may consider the problem of description of Finsler metrics 
whose geodesies are circles of a fixed geodesic radius R on the unit sphere. 
TheoremlHland Corollary El still apply, see Remark 13.31 However the situation 
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is different on S"^, as far as the continuous functions are concerned whose in- 
tegrals vanish over all geodesic discs of radius R. For all but countably many 
special values of R, such functions are identically zero, see ^29j. This implies 
an interesting "almost everywhere" rigidity: for a generic R, there is only 
one (up to summation with exact 1-forms) Finsler metric whose geodesies 
are circles of radius R; this unique metric is an analog of the metric (jH} in 
R^. Of course, in the plane, all values of the radius R are equivalent, due to 
similarity. 
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